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Using 14M ip(2S) events accumulated by BESII at the BEPC, a Covariant Helicity Amplitude Analysis is 
performed for tp(2S) — > tt + tt~ J/ip, J/ip — > n + . The 7r + 7r~ mass spectrum, distinctly different from phase 
space, suggests a production in this process. Two different theoretical schemes are used in the global fit to 
the data. The results are consistent with the existence of the a. The a pole position is determined to be 



1. Introduction 

We report here a study of the process tp(2S) — ► 
7r + 7r~ J/tp, which is the tp(2S) decay mode with 
the largest branching fraction [1], using very clean 
tp(2S) — > TT + ir~J/if.> (J/tp — > events. 
Early investigations of this decay by Mark I [2] 
found that the tt + tt~ mass distribution is strongly 
peaked toward higher mass in contrast to what is 
expected from phase space. Furthermore, angular 
distributions favored S-wave production between 
the tttt system and J/tp, as well as an S-wave de- 
cay of the dipion system. 

BESI studied this process with much higher 
statistics (3.8 million tp(2S) events) and found 
that an additional small D-wave component was 
required in the decay of the dipion system [3]. 
Also various heavy quarkonium models were fit- 
ted, and the parameters for these models deter- 
mined [3]. 

Here, we fit the tt + tt~ system from tp(2S) — * 
ir + Tr~J/ip decays with the J — ++ a me- 
son. In this decay, the interaction between the 
tttt system and ip(2S) or J/tp is small since these 
charmonium states are very narrow, so the dipion 
system is a quasi-isolated system [4] . 

The a meson was introduced theoretically in 
the linear a model [5], and its existence was 
first suggested in a one-boson-exchange potential 
model of nuclear forces [6]. The a meson is im- 
portant due to its relation with dynamical chiral 
symmetry breaking of QCD [7]. 

There was evidence for a low mass pole in early 
DM2 [8] and BESI [9] data on J/tp -> lotttt. A 
huge event concentration in the 7 = S-wave 
tttt channel was observed in a pp central produc- 
tion experiment in the region from m wn = 500 
to 600 McV/c 2 [10]. This peak is too large to 
be explained as background [11]. Many stud- 
ies on the possible resonance structure in tttt 
scattering have appeared in the literature [12]. 
It was proved that the existence of a light and 
broad resonance is unavoidable even with non- 
linear realization of chiral symmetry [13]. Care- 
ful theoretical analyses were made to determine 
the pole location, which was found to be M — 
iT/2 = (470 ± 30) - i(295 ± 20) MeV/c 2 [14] and 
M-iT/2= (470±50)-i(285±25)MeV/c 2 [15]. 



Renewed experimental interest arose from E791 
data on D + — ► tt + tt~tt + [16], where it was found 
that M = 478±23 ± 17 MeV/c 2 , T = 324^ ± 
21 MeV/c 2 . In the recent partial wave anal- 
ysis of the decay J/tp — > lott + tt~ [17] by BE- 
SII, the pole position of a was determined to be 
(541 ± 39) - i(252 ± 42) McV/c 2 . All these exper- 
imental results still have large uncertainties. 

Fig. 1 shows the decay mechanism of tp(2S) — > 
7r + 7r~ J/tp in the S-matrix formalism. There are 
three main contributions including an S'-wave res- 
onance (a), a D-w&ve term (2 + ), and a con- 
tact term which is the destructive background re- 
quired by chiral symmetry [18]. The total am- 
plitude is the sum of these three components. 
The decay tp(2S) — > TT + n~J/tp can also be de- 
scribed with an effective Lagrangian for the vector 
pseudo-scalar pseudo-scalar (VPP) vertex, along 
with the tttt final state interaction (FSI) obtained 
from tttt scattering data in a Chiral Unitary Ap- 
proach (ChUA) [19]. In such an approach, the 
a resonance is generated dynamically as a pole 
of the unitarized i-matrix, and the pole position 
is 469 - i203 MeV/c 2 [20]. A fit to the BESI 
tp(2S) — ► ir + 7r~J/tp data shows that the tttt FSI 
plays an important role in this process [19]. A 
similar result was obtained in Ref. [21] with a 
comparison of the cases with and without the tttt 
FSI. We fit our data with both the S'-matrix and 
ChUA schemes. 

2. BESII Experiment 

The data sample used for this analysis is taken 
with the BESII detector at the BEPC storage ring 
operating at the tp(2S) resonance. The number of 
tp{2S) events is 14.0±0.6 million [22], determined 
from the number of inclusive hadrons. 

The Beijing Spectrometer (BES) detector is 
a conventional solenoidal magnet detector that 
is described in detail in Ref. [23]; BESII is the 
upgraded version of the BES detector [24]. A 
12-laycr vertex chamber (VC) surrounding the 
beam pipe provides track and trigger informa- 
tion. A 40-layer main drift chamber (MDC), lo- 
cated radially outside the VC, provides trajectory 
and energy loss (dE/dx) information for charged 
tracks over 85% of the total solid angle. The mo- 
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Figure 1. Decay mechanisms for ip(2S) — ► tt + tt J/i/j in the S'-matrix formalism. The final amplitude is 
the superposition of ip(2S) — > <r <//?/>, ip(2S) — > 2 + J/tp, and ip(2S) — > (7'' + 7'' _ )cont'-V , /' ^-matrix elements. 



mentum resolution is <7 C /p = 0.017^1 +p 2 (p 
in GeV/c), and the dE/dx resolution is ~ 8%. 
An array of 48 scintillation counters surrounding 
the MDC measures the time-of-flight (TOF) of 
charged tracks with a resolution of ~ 200 ps for 
hadrons. Radially outside the TOF system is a 12 
radiation length, lead-gas barrel shower counter 
(BSC). This measures the energies of electrons 
and photons over ~ 80% of the total solid angle 
with an energy resolution of <je/E = 22%/\[E (E 
in GeV). Outside the solenoidal coil, which pro- 
vides a 0.4 Tesla magnetic field over the tracking 
volume, is an iron flux return that is instrumented 
with three double layers of counters (MUID) that 
identify muons of momentum greater than 0.5 
GeV/c. 

A GEANT3 based Monte Carlo (MC) simu- 
lation program [25] with detailed consideration 
of detector performance (such as dead electronic 
channels) is used to simulate the BESII detector. 
The consistency between data and MC simula- 
tion has been carefully checked in many high pu- 
rity physics channels, and the agreement is quite 
reasonable [25]. 

3. Event Selection 

Events with 7r + 7r~/i + p~ final states and with 
the invariant mass m^+^- constrained to the J /ip 
mass are selected for analysis. Each track, recon- 
structed using hits in the MDC, must have a good 
helix fit in order to ensure a correct error matrix 



in the kinematic fit, and the number of tracks are 
required to be between 4 and 7. 

To select a pair of muons, the muon pair can- 
didates tracks are required to have net charge 
zero; p^+ > 1.3 GeV/c or p^- > 1.3 GeV/c or 
2fy+ + Pp.- > 2.4 GeV/c; |cos0 M | < 0.6 to en- 
sure that tracks are in the sensitive region of the 
MUID; the cosine of the angle between these two 

tracks in their rest frame cos^f 1 !? _ < —0.975 to 

fj, fi 

guarantee the collinearity of the tracks; the sum 
of the MUID hits N% u + N Mt > 3 to ensure that 
the tracks are muons; and the invariant mass of 
two candidate tracks rrin+u- within 0.35 GeV/c 2 
of the J /tp mass. 

For 7r + 7r~ pair selection, the two candidate 
tracks are also required to have net charge zero. 
Each track is required to have momentum p n < 
0.5 GeV/c, polar angle | cos#| < 0.75, and trans- 
verse momentum p nxy > 0.1 GeV/c to reject 
tracks that spiral in the MDC. The dE/dx mea- 
surement of each track must be within three stan- 
dard deviations of the dE/dx expected for the 
pion hypothesis, and the cosine of the laboratory 
angle between the candidate tracks must satisfy 
cos^tttt < 0.9 to eliminate e + e~ pairs from 7 con- 
versions. The mass recoiling against the candi- 
date 7r + 7r~ pair, m r ^+™L, is shown in Fig. 2a. 
In order to get well reconstructed signal events 
and to suppress background, \m r ^° %1 : — mj/^\ < 
20 MeV/c 2 , corresponding to three times the 
mass resolution, is required. 
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With the above selection criteria, about 40,000 
i>{2S) — > 7r + 7r~J/-0 — > 7r + Tr~fi + (i~ candidate 
events are obtained. Fig. 2b shows the 7r + 7r~ in- 
variant mass distribution for these events, where 
the dots with error bars are data, and the his- 
togram is Monte Carlo simulation with the PP- 
GEN generator, which is based on chiral symme- 
try arguments and partially conserved axial vec- 
tor currents [26]. It describes the low mass tttt 
spectrum reasonably well but not the high mass 
region; the inconsistency between data and Monte 
Carlo will be considered in the systematic errors. 

The main background channels are 
^(25) -> r,J/i> (J/4> -> /x+M"), 1>(2S) -> 
ir+Tr-J/4* (J/ip -> TT+TT-) and tp(2S) -> 
n + ir~J/ip {J/ip — * /twr)- However, MC simulation 
indicates that their total contribution is less than 
0.1%, which can be neglected. The contamination 
from continuum production e + e~ — > tt + tx~ ji + [i~ 
is also very small and neglected in this analysis. 

4. Analysis Method 

Two different schemes are used to fit our data. 
In the first, based on the diagrams in Fig. 1 and 
taking the VPP vertex as a constant, the total 
differential cross section which describes ij>(2S) — > 
tt + tx~ J/ip is 



da 



MA„ 



E 

MX* 



A + E ^ 2 + A 



contact 



(1) 



where A s represents the amplitude for ip(2S) — > 
XJ/ip — > Tr + ir~J/ip with the spin of X being s, 

represents the solid angle, M is the magnetic 
quantum number along Z-axis of iji(2S), and 
A2 are the helicities of the J ftp and 2 + compo- 
nents, respectively, and A confa ct is the amplitude 
of the contact term. 

In the second, considering the VPP vertex and 
the S'-wave tttt FSI, while neglecting the D-wave 
FSI, the amplitude is [19] 



A = Vq+ Vqs G-2t 



1=0 

"7T7T — >7Z7T ) 



(2) 



V = --72 (fflPl ' P2 + ff2P?P° + 53«4) e * ' 6 > ( 3 ) 



where G is the two-pion loop propagator, Vos 
is the S-wave part of Vo, and t^^ n is the full 
S-wave 1 = tttt — > 7T7t i-matrix, which is the 
same as those defined in Refs. [19,20]; p\ and P2 
are the four momenta, of the two pions, and p\ 
and P2 are their energies in the lab frame; c/i, (72, 
and 33 are free parameters to be determined by 
data. 

The normalized probability density function 
used to describe the whole decay process is 



f(x,a) = 



da/dfl 



(4) 



where x represents a set of quantities which are 
measured by experiment, and a represents un- 
known parameters to be determined. The total 
cross section, a, can be expressed as 



(5) 



where e(£l) is the detection efficiency which is usu- 
ally a function of detector performance. The total 
cross section can be determined by MC integra- 
tion. Re-weighting a total of TV generated events 
based on simulated ip(2S) — > 7r + 7r~ J/ip using a 
phase space generator, the total cross section is 
then 
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N„ 



E c J da 
I WO 



(6) 



where N mc (< N) is the number of MC simu- 
lated events after applying the selection criteria. 

The maximum likelihood function [27,28] is 
given by the joint probability density of the se- 
lected ip(2S) — ► 7r + 7r~/! + /i~ events, 



N evt 

£=H f(x,a) 



(7) 



and a set of values, a, is obtained by minimizing 
the function 5, 



S 



log£ 



(8) 



For the amplitudes in the first model, the am- 
plitudes for the cascade two-body decay process 
can be expanded with helicity amplitudes as: 



A s = FLD 



. v BW x (Sn,m x ,Tx)F5oD 



A,0 ) 



(9) 
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Figure 2. Distributions of candidate ip(2S) — ► 7r + 7r~ J/ip (J/ip — > /i + /U~) events, (a) is the 7r + 7r~ recoil 
mass spectrum, fitted with a double Gaussian function, and (b) is the 7r + 7r~ invariant mass spectrum. 
The histogram in (a) is data, the curve is the fit, the events between arrows are selected; in (b) dots with 
error bars are data, and the histogram is MC simulation. 



where F^ v is the helicity amplitude, which can 
be found in Ref. [28], D J MX _ v (<j) t 6,Q) is the D- 
function, and BWx{S nv , n%x , Tx) is the Breit- 
Wigner propagator of X, defined as : 



BW x (s^,m x ,Tx) 



- m\ + im x T x {s^) 



(10) 



The a particle, a broad structure in the low 
7r + 7r~ mass region, is not a typical Breit-Wigner 
resonance. In the first model, four types of Breit- 
Wigner parameterizations are used to describe it: 

1) Constant width 



3) P.K.U. ansatz [29], which removes the 
spurious singularity hidden in Eq. (12) 



r *( s ) = PZjT T = \ 1 - 



•x 



Ami 



(13) 



'A' 



4) Zou and Bugg's approach [30], where 
the form includes explicitly into Tx (s) the 
Adler zero at s = m\/2. 



fji 



7 2~T + 92 7 2~T ( , i4 / 



s-ml/2 
9i = f{sj— — e 



ml/2 



T x (s) = r. 



(11) 



where the definitions of p X7r , pi n , and /(s) 
are the same as Ref. [17]. 



2) Width containing a kinematic factor, 
which was used by the E791 Collabora- 
tion [16] 



T x (s)=pT=\l 



4m| 



(12) 



5. Partial wave analysis 

The minimization used in the partial wave 
analysis and to obtain the pole parameters of the 
a is based on MINUIT [31]. For the first model, 
the components considered include amplitudes of 
ct(0 + ), a D-wave term, and a contact term. The 



G 



tail of the / (980) has been tried in the fit. How- 
ever, it has similar behavior to the contact term 
in this mass region, and therefore it is ascribed 
to the contact term. All four a Brcit-Wigner pa- 
ramctcrizations fit the data well, but have strong 
destructive interference with the contact term, es- 
pecially in the low 7r + 7r~ invariant mass region. 
The D-wave contribution is only 0.3 to 1%, in 
agreement with the BESI result [3] based on a 
different analysis method. Fig. 3 shows the pro- 
jections of the fit results compared with data for 
the Breit-Wigner parameterization for the P.K.U. 
ansatz; other parameterizations give similar re- 
sults. 

The global fits determine the best estimation of 
the Brcit-Wigner parameters for each parameter- 
ization. The pole position in the complex energy 
plane is related to the mass and width of the res- 
onance by 

^Spole = m a - i-^- . (15) 

The best fit results and the corresponding pole 
positions for all the parameterizations are listed 
in Table 1. The statistical error of the resonance 
mass (width) is determined by a decrease of \ in 
the log-likelihood from its maximum value with 
all other parameters fixed to their best solutions. 

For the second model, where the VPP vertex 
is represented by an effective Lagrangian and the 
7T7T 5*- wave FSI is included, the 7T7t mass spectrum 
can also be reproduced well. Here the a requires 
a much smaller interference between the S-wave 
FSI and the contact term. In this case, the fit is 
worse than the fits of the first model; this may 
due to the fixed pole position of the a and the 
neglected D-w&ve contribution. For the second 
model, the pole is not measured in the fit, but 
taken from Ref . [20] , which was determined from 
7T7T scattering data. 

To check the goodness of fit in our analysis, we 
construct a variable 

where N is the number of cells, Np T and 

are the numbers of events in the i'th cell of the 



Dalitz plot with axes m 2 , _ and m 2 , , , for data 
and MC simulation, respectively. Such a vari- 
able should be distributed according to the \ 2 
distribution with n = N — K degrees of freedom, 
where K = 12 is the number of parameters to 
be determined in our Maximum Likelihood fit. 
In our case, 15 bins in both m 2 + _ and m 2 T/ , , 

give 225 cells. To ensure proper x 2 ( n ) behavior, 
cells with less than five events have been merged 
into adjacent ones. The number of cells becomes 
N = 208, and the number of degrees of freedom 
n = 196. From the observed \ 2 value determined 
using Eq. (16) for each parameterization, the con- 
fidence levels (C.L.) arc calculated and listed in 
Table 1. 

6. Systematic Errors 

For the first model, the systematic error of the 
a pole position arises from the uncertainties of 
the strength of the 2 + component, the form of 
the contact term, and the inconsistency between 
data and of MC simulation. For the 2 + compo- 
nent, we conservatively remove it from the fit, and 
the difference of the fitted values from the nom- 
inal values are taken as systematic errors. Two 
contact terms, namely, constant amplitude and 
a.\ + ict2P, where a\ and «2 are two parameters 
to be fitted, are adopted in the fit; the difference is 
considered as the systematic error. The MC sim- 
ulation and data have different mass resolutions 
in the high mass region of the 7r + 7r~ system. A 
modification of 7r + 7r~ mass resolution is made to 
improve the fit, and the difference of the fitted 
pole positions with and without this modification 
is taken as the systematic error. The systematic 
error from non-signal backgrounds is neglected. 

In order to obtain the m CT and -# errors in 
Eq. (15), we set the denominator in Eq. (10) 
equal to zero and obtain the pole position and 
corresponding errors by taking into account the 
mass and width errors of the Brcit-Wigner pa- 
rameterizations. This is done using a MC sam- 
pling method, where the correlation between the 
mass and width is ignored. Table 2 summarizes 
the systematic errors from all sources, and Table 
1 lists the parameters of pole position and their 
total errors. 
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Figure 3. Fit results of ^(2S I ) — > ir + ir~ J /if) (P.K.U. ansatz). Dots with error bars are data and the 
histograms are the fit results, (a) and (b) are the ir + Tt~ invariant mass, (c) and (d) the cosine of the a 
polar angle in the lab frame, and (e) and (f) the cosine of the 7r + polar angle in the a rest frame. The 
upper plots are the detected distributions, while the bottom ones are the distributions after efficiency 
correction. 




Figure 4. The 7r + 7r~ invariant mass distribution including the components. The left is the P.K.U. ansatz, 
which has the contributions from a, D-wave term, contact term, and their sum (the Z?-wave is enlarged by 
a factor of 20 in the figure). The right is the tt + tt~ invariant mass fitted by the formula from Ref. [19,20], 
with no explicit D-wave. Dots with error bars are data, and the histograms are the fit results. 
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Tabic 1 



Fit results of the two models and all the Brcit-Wigner parameterizations 



Model 


Constant T 


r with p 


P.K.U.ansatz 


Bugg &; Zou's approach 


Ref. [19,20] 




(553 ± 15 ± 47) 


(559 ± 6 ±26) 


(554 ± 13 ± 65) 


(541 ± 9 ± 95) 


469 - i203 


Pole(MeV/c 2 ) 


-j(254 ± 23 ± 54) 


-i(179 ± 7 ±18) 


-£(240 ± 4 ± 19) 


-i(253 ±8±33) 


(input) 




140308 


72735 


133208 


171586 


30765 


-^contact t 


121625 


63133 


111230 


157741 


3039 


-lnL 


-16174.5 


-16166.8 


-16171.0 


-16174.3 


-15974.1 


xl bs Jndf 


217.83/196 


227.54/196 


224.07/196 


217.88/196 


392.73/208 


C.L. 


0.1362 


0.0608 


0.0825 


0.1357 


3 x 10~ 13 



f : N a and N contact are the numbers of events in the fit. 



Table 2 

Systematic errors in the pole position ( MeV/c 2 ) 





Constant F 


r with p 


P.K.U. 


ansatz 


Bugg & 


; Zou's approach 






T/2 


Cm 


or/2 




o"r/2 


Cm 


o"r/2 


2 + uncertainty 


22 


35 


4 


14 


10 


2 


3 


3 


form of contact term 


11 


14 


1 


2 


30 


8 


91 


22 


M.C. imperfection 


40 


38 


25 


10 


56 


17 


37 


24 


Total Error 


47 


54 


26 


18 


65 


19 


95 


33 



7. Results and discussion 

The process ip(2S) — > iy + tt~ J/ip, J/ip — > 
/i + yT is studied based on 14 x 10 6 ip(2S) events 
collected with the BES11 detector. The 7r + n~ 
invariant mass spectrum of ip{2S) — ► 7r + 7r _ J /if) 
has a severe suppression near the tt + tt^ thresh- 
old, which is distinctly different from phase space 
and suggests a production in the process. We fit 
the data with two different models. For the first 
model, using four different Breit-Wigner parame- 
terizations, the data can be well fitted, although a 
strong cancellation between the a and the contact 
term is required. In fact, such a large cancellation 
is dictated by chiral symmetry [26,18]. The pole 
positions of a are determined for different Brcit- 
Wigner parameterizations, which are (553 ± 15 ± 
47) - i(254 ± 23 ± 54) MeV/c 2 (constant width), 
(559 ± 6 ± 26) - i(179 ± 7 ± 18) MeV/c 2 (width 
containing kinematic factor), (554 ± 13 ± 65) — 
i(240 ± 4 ± 19) MeV/c 2 (P.K.U. ansatz), and 
(541 ± 9 ± 95) - i(240 ± 8 ± 33) MeV/c 2 (Bugg 
& Zou's approach). The first Breit-Wigner pa- 
rameterization may be problematic because the 
imaginary part does not vanish at threshold. The 
second parameterization gives a small a width, 
and creates a virtual state in the real energy axis 
below the tttt threshold [29]. The final best esti- 



mate of the a pole position from this analysis is 
(5521^) - i(232±^) MeV/c 2 , where the central 
values are obtained by a simple mean of the differ- 
ent Brcit-Wigner parameterizations excluding the 
second one, while the errors cover the statistical 
and systematic errors, including the differences in 
the three Breit-Wigner parameterizations. 

We also fit our data according to the scheme in 
Ref. [19]. It is found that the 7T7t S-wave FSI plays 
a dominant role in '0(2S') — > tt + tt~ J j 'ijj , while the 
contribution from the contact term is small. This 
means that the a meson has a significant contri- 
bution in this process. The a pole used in this 
fit, 469 - i203 MeV/c 2 is consistent with the fits 
to the Breit-Wigner functions. This implies that, 
although the two theoretical schemes are very dif- 
ferent, both of them find the a meson at similar 
pole positions. 

If the a meson exists, the pole should occur 
universally in all 7T7t system with correct quan- 
tum numbers. Our analysis demonstrates that, 
in ^(25) — > 7T + 7r- J/V^, even though there is no 
apparent peak structure, one can still determine 
the pole location in good agreement with that 
obtained from J/ip — > lutt + tt^ decay [17] by as- 
suming a simple form of the contact term. Hence 
it provides further evidence for the a meson. 
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